Flow effects on the freeze-out phase-space density in heavy-ion collisions 
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The strong longitudinal expansion of the reaction zone formed in relativistic heavy-ion collisions 
is found to significantly reduce the spatially averaged pion phase-space density, compared to naive 
estimates based on thermal distributions. This has important implications for data interpretation 
and leads to larger values for the extracted pion chemical potential at kinetic freeze-out. 
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The phase-space density of mesons produced in ultra- 
relativistic heavy ion collisions is an interesting quantity. 
Its spatial average at freeze-out (symbolised by t{), 



(f)(p) 



J d 3 x f 2 (t>t f ,x,p) 
J d 3 xf(t>t f ,x,p) 



(1) 



can be measured by combining the single-particle mo- 
mentum spectrum with the "homogeneity volume" ex- 
tracted from Bose-Einstein correlation measurements [Q . 
Bertsch's original formula El was refined in S to ex- 
clude contributions from long-lived resonances decaying 
far outside the collision fireball, and a relativistically co- 
variant derivation was given in |],^| . With these improve- 
ments Bertsch's formula reads 



</>(p) 



y/Mpj dn/ (dyp±dp± d<p) 



(E p /n^)R s (p) ^Rl{p)$i{p)-RUp) 



(2) 



Here A is the intercept at vanishing relative momentum 
of the Bose-Einstein correlation function, and the factor 
y/ X(p) corrects for contributions from decays of long- 
lived resonances in accord with the core-halo model 
The numerator is the (thus corrected) invariant momen- 
tum spectrum, and the denominator contains the ho- 
mogeneity volume calculated from the Hanbury Brown- 
Twiss (HBT) radii extracted from the Bose-Einstein cor- 
relation function Q). 

The strongest motivation for measuring the phase- 
space density of pions at freeze-out (i.e. on the last- 
scattering surface) comes from the search for multi-boson 
symmetrization effects Q. Analyses of the hadron yields, 
spectra, and two-particle correlations from Pb+Pb col- 
lisions at the CERN SPS @ and Au+Au collisions at 
RHIC H indicate that freeze-out occurs in two stages: 
chemical freeze-out, where the particle yields decouple, 
happens first, reflecting a temperature of around 170 
MeV, while kinetic or thermal freeze-out takes place 
much later, at temperatures around 100-120 MeV. Be- 
tween these two points the system cools adiabatically at 
constant particle numbers which requires the build-up 
of positive chemical potentials M. It has been specu- 
lated that for pions, with only 139 MeV rest mass, this 



chemical potential could approach the Bose condensa- 
tion limit ^ 7r =m 7r . If this were true, measurements of 
the above type should exhibit large values of (f)(p) at 
small p. Large phase-space densities also generate sig- 
nificant multi-boson symmetrization effects on the pion 
spectra 10 and Bose-Einstein correlation functions (see 
|TTf and references therein), which lead to a reduction 
of the homogeneity volume extracted from the standard 
HBT correlation analysis Q] . This should cause an addi- 
tional enhancement of (f)(p) as determined by Bertsch's 
formula (^J), rendering it at the same time unreliable as 
an estimator of the real average phase-space density. 

Previous analyses of heavy-ion data at the AGS j^j 
and SPS |l~2|,|3| concluded that there was no evidence for 
unusually large pion freeze-out phase-space densities at 
these collision energies. This conclusion was based on a 
comparison of (/) (p) from Eq. (g) with a Bose-Einstein 
distribution, yielding rough agreement when inserting a 
kinetic freeze-out temperature of around 100-120 MeV 
(as extracted by other methods ) and a small or vanish- 
ing pion chemical potential. In some cases [|3|]T3]] compar- 
ison was made with a transversely boosted Bose-Einstein 
distribution, in order to account for transverse flow ef- 
fects. We here point out that all these comparisons may 
have been misleading since they neglected a strong reduc- 
tion effect on the spatial averaging originating from the 
longitudinal expansion of the source. We show that for 
longitudinally expanding sources the spatially averaged 
phase-space density is smaller than its thermal value in 
the local rest frame, and that accounting for this effect 
is likely to considerably increase the extracted value of 
the pion chemical potential at kinetic freeze-out. In this 
paper we concentrate on a theoretical exposition of the 
basic mechanism and leave a (re) analysis of existing data 
for later. 

For a qualitative understanding of the effect under dis- 
cussion, let us for simplicity assume a fireball of constant 
matter density. We fix the momentum and focus on the 
averaging over position space. If the (thermalized) fire- 
ball expands, different parts moving with different ve- 
locities contribute to the production of particles with 
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fixed momentum p with different rates, given by Bose- 
Einstein distributions boosted by the local flow velocity. 
The phase-space density at momentum p then becomes 
a function of position space. For an average of type (|l]) , 
where the function is weighted with itself, the result is 
always smaller than the maximum of the function (here 
corresponding to zero flow velocity). In general the p- 
dependence of the spatial average will no longer be of 
Bose-Einstein form. [Obviously, there is no position de- 
pendence of the emission rate in a non-expanding source 
for which the spatial average would thus coincide with 
the underlying distribution.] 

We illustrate this effect quantitatively within a spe- 
cific model for the expanding source. We will assume a 
thermalized and longitudinally boost-invariant fireball. 
This simple model is expected to provide a good de- 
scription of pions emitted near midrapidity at RHIC 
and LHC energies. We show that the calculated ^in- 
dependence of the spatially averaged phase-space density 
is flatter than a Bose-Einstein distribution with the as- 
sumed freeze-out temperature. The major contribution 
to this flattening arises from the longitudinal expansion 
which strongly reduces (/) (p) relative to the correspond- 
ing static Bose-Einstein distribution. Transverse flow in- 
troduces a weaker additional flattening by shifting weight 
from low to high transverse momenta p±. 

The phase-space density f(t,x,p) is related to the 
source emission function S(x,p) by 



f(t,x,p) = 



En 



dt' S{t',x + v(t' -t),p), (3) 



where v =p/E p is the velocity corresponding to the mo- 
mentum p. The prefactor (2ir) 3 /E p ensures the standard 
normalisation 



N = 



W /ft *' p) = /f d ' xS{x ' p) - (4) 



We take N as the average number of "directly emitted" 
pions, including those from the decays of shortlived res- 
onances but excluding those from resonances with l ong 
lifetimes For simplicity, we follow Refs. [0JT| 

and assume that the change in shape and normalisation 
of the thermal pion distribution at freeze-out caused by 
adding pions from shortlived resonance decays can be ef- 
fectively absorbed by giving the pions a positive chemical 
potential. Although this cannot replace a full resonance 
decay calculation using the proper decay kinematics [^6) , 
it qualitatively reproduces the decay-induced low-pj_- 
cnhancement of the measured pion spectra 14 16 j. 

Eq. (0) is brought into the form (^) by following the 
steps outlined in Sec. 3.4 of Ref. Q. In order to include 
all "directly emitted" pions, the time t on the right hand 
side of Eq. (Q) must be later than the time at which the 
last pion was emitted. Due to Liouville's theorem, the 
average phase-space density ([l]) remains constant after 



completion of the freeze-out process. This allows [|| to 
replace in ([!]) the spatial integrals at constant time t by 
integrals over the freeze-out hypersurface £f(x) describ- 
ing the last pion scattering or production points: 



</>(p) 



J Sf p ■ d 3 a(x) f 2 (x,p) 
/ E( p- d 3 a(x) f(x,p) 



(5) 



If Sf(x) is parametrised by an cc-dependent freeze-out 
time tf(x), then 



p-d*o-(x)r(t,x, P ) = 

E p J d 3 x(l-vV x t { (x)) f n (t f (x),x,p). (6) 

For the model emission function we take jl7) 
m± cosh(y— rj) S(t—T{) 



S( X >P) ,o_ys | i> u{x)—fi(r) "1 j 



(2k) 



(7) 



exp 



where 



p-u(x)=m± cosh(y—rj) cosh^(r) — p± cos(4>—ip) sinh£(r). 



We use longitudinal proper time T=\/t 2 —z 2 , space-time 
rapidity rj = | ln[(t+z)/(t— z)], and transverse polar co- 
ordinates (r, ip) to parametrise x. The pion momentum 
is p = (m± cosh y, p± cos tp, p± sin tp, m± sinh y) . Eq. (Q) 
describes a longitudinally infinite source with boost- 
invariant longitudinal expansion and sharp freeze-out at 
proper time T{ . Transverse expansion is parameterised by 
the transverse flow rapidity C( r )i to be specified below. 
The r-dependent chemical potential fj,(r) allows to dis- 
cuss different transverse density profiles of the source at 
freeze-out. While probably too simple for a quantitative 
comparison with data, this model is technically conve- 
nient and allows to investigate the effects of longitudinal 
and transverse flow as well as the influence of the chem- 
ical potential on the average phase-space density (f)(p). 

We now evaluate Eq. (@) for midrapidity pions (y = 0). 
We can rotate the coordinate frame such that (j> = 0. 
From (^) we obtain for the phase-space density imme- 
diately after freeze-out 



f(t{(x),x,p) 



exp 



p ■ u(xf)-n(r) 
T 



(8) 



with Xf = (tf(x), x). For the integration in Eq. we use 
the freeze-out surface at r = T£ with integration measure 

(l — v-V x tf(x)) d 3 x = Tf cosh rj dr\ rdrd<p. (9) 

To obtain (f)(p) from (||) we need (E p — mj_ for y = 0) 

cosh rjdrjr dr dip 



A n {p±) = m±Ti 



exp 



/ p-a(gf)— n(r) \ 
\ T J 



1 



(10) 
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for n = l,2. Note that Ai(p±) is, up to a factor (2ir)~ 3 , 
the invariant momentum spectrum E p (dN/d 3 p). 

We evaluate A n by expanding the Bose distribution as 
(e x — = J2T=i e ~ kx - This reduces ( |l0| ) to a sum of 
exponential integrals. We now integrate over 77 and ip in 
the standard way JL8| , obtaining Bessel functions Ki and 
Io, respectively. Only the radial integration must be done 
numerically With the shorthands 



aj_(r) = p± sinh C(r)/T 
x (r) = m ± coshC(r)/T 



we find 



(/) (Px) 
Mpx) 

A 2 (p±) 



fc=i 



x;(*-i)^(p±), 

/■CO 

47TTO ± r f / rdre fc[Al(r)/TI 
./o 

xI (feax(r))K 1 (fc ( 5 ± (r)) . 



(11a) 
(lib) 



(12a) 
(12b) 

(12c) 
(12d) 



We will consider two models for the transverse density 
profile: a box-like density 



[box] fi(r) = 



[i B for r < i? box 
—00 otherwise 



and a Gaussian profile 



[Gauss] /i(r) = fiQ — T 



(13a) 



(13b) 



Here \ib and \iq are the chemical potential values in the 
center of the fireball. For the Gaussian profile, fie will 
be larger than the global (averaged) chemical potential 
used in chemical analyses of particle abundances p9| . 

As a warm-up, let us consider the case without trans- 
verse flow, £(r) = 0, in the Boltzmann approximation (i.e. 
keeping only the first term in the sums over k). Then 



</?=T) (Px) = ^ 



where A,r is the spatially averaged pion fugacity 



J rdre 2»(r)/T 

Jrdr e^/ T ' 



(14) 



(15) 



the average being taken with the transverse matter den- 
sity p(r) ~ e^ r 'l T . A box-like density gives A^ : 



whereas a Gaussian profile leads to A^ : 



In Fig. we show the result (JLJ) for T= 120 MeV and 
three values for the average pion fugacity A^ . According 
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FIG. 1. Average phase-space density (/) (p±) at rapidity 
y = for pions from a longitudinally expanding fireball with 
Boltzmann-distributed particle momenta, for a temperature 
of 120 MeV and three representative values for the pion chem- 
ical potential. 



to (|lj) the curves do not depend on the transverse den- 
sity profile. The calculation in Fig. [j] was done with a 
box-like density profile for which fis can be considered 
as the global chemical potential. The same curves would 
be obtained in the Gaussian model by adjusting fiQ such 
that e MG / T = 2 e^ B l T (leading to a twice larger fugacity 
X n [r—0) — 2\ lx in the fireball center). 

The dotted line in Fig. |l] is shown for reference and 
represents a Boltzmann distribution with vanishing pion 
chemical potential for the same temperature as used in 
the calculations (i.e. the original thermal pion momen- 
tum distribution in the local rest frame). Comparing it 
with the solid curve one sees that, even without trans- 
verse flow, at low p± the spatially averaged phase-space 
density of the longitudinally expanding source is about 
35% lower than its value in the local rest frame. In order 
to obtain an (/) (p±) which looks like a Boltzmann distri- 
bution with vanishing pion chemical potential, we need 
in fact an average pion fugacity A w « 1.55. For practical 
purposes this interpretation is better presented the other 
way around: if the curve corresponding to A^ w 1.55 (or 
Hb = 50 MeV at T= 120 MeV) was measured and then 
fitted by a Boltzmann distribution with T= 120 MeV, 
one would erroneously conclude /is — 0. 

For large m x , Ki(2m_ L /T)/Ki(m_ L /T) » e - m ^/ T , and 
this type of error no longer arises. Unfortunately, this re- 
quires for pions much larger values of px than presently 
experimentally accessible. The reason for the strong 
reduction at low p±, resulting from taking the spatial 
average, is that this average extends over a large lon- 
gitudinal homogeneity length, and that for a longitu- 
dinally expanding source f(x;y=0,p±) agrees with the 
static Boltzmann distribution exp(— m±/T) only at the 
point r\ = 0, being smaller everywhere else by a factor 
exp[— m^(cosh n — 1)]. For large p± the longitudinal ho- 
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mogeneity length decreases [H and the reduction due to 
spatial averaging becomes less severe. 
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FIG. 2. Average phase-space density (/) for pions with 
y — from a longitudinally and transversally expanding fire- 
ball with Bose-Einstein distributed momenta of temperature 
120 MeV. The transverse flow is controlled by rj t . The trans- 
verse density profile has the shape of a box. 

We now discuss the more general case of a source 
with additional transverse flow, including the full Bose- 



Einstein distribution by summing over all k in Eqs. (12) 



The transverse flow profile £(r) influences the radial av- 
eraging of the factors containing the pion chemical po- 
tential, and it is no longer possible as in ([u]) to factorize 
an effective average pion fugacity A w from the result for 
(/) (P-i-)- For the transverse flow rapidity we take a linear 
profile with slope rjt, 



where r rms = R hox /V2 and r rms = \/2 i?Ga USS , respective- 
ly, are the rms radii corresponding to the two different 
transverse density profiles. Thus r\ t = C rms can be inter- 
preted as the average (rms) transverse flow rapidity. 

Figures]^ and|| show that the longitudinal flow is again 
responsible for most of the suppression of (f)(p) com- 
pared to the local rest frame distribution. Transverse 
flow tends to flatten the pj_-dependence of (/) rather than 
to suppress (/). This effect stems from enhanced con- 
tributions at higher transverse momenta from outward- 
moving fluid cells; this is the well-known "blueshift" of 
the transverse momentum spectrum. The details of this 
flattening depend on the particular model for the trans- 
verse density profile. 
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FIG. 3. Same as Fig. ^ but for a Gaussian transverse den- 
sity profile. 

The calculated (f)(p) in Figs. || and || are compared 
to Bose-Einstein distributions with an averaged chemical 
potential /2 which takes into account the averaging over 
the transverse density profile. It is obtained by calculat- 
ing the "transverse average" 



/tr(p) = 



C(r) = Vt 



(16) 



f™rdrf tI (r,p) 

of the underlying local distribution without flow 

1 

ftr(r,p) 



exp [(m± - fJ.(r))/T] - 1 



(17) 



(18) 
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where (i(r) characterizes the chosen transverse density 
distribution. Then a Bose-Einstein distribution /be = 
l/[exp((mj_ — p)/T) — 1] with the same temperature is 
fitted to /tr(p), using p as a fit parameter to ensure the 
same normalisation J fBE{fn±)n r i'±d m _L." Obviously, for 
the box profile p — [Xb'-, for the Gaussian profile, p < /ic- 
For the Gaussian profile the Bose fit is reasonable but 
not perfect since ftr(p) is a bit more strongly peaked 
at small p± than for the box profile. However, since 
p/T controls the normalization of /, i.e. the total pion 
yield at midrapidity y — 0, p is identical with the global 
chemical potential extracted from a chemical analysis of 
midrapidity particle abundances at the thermal freeze- 
out temperature, and therefore a good reference value. 

For the direct comparison of the two transverse den- 
sity profiles we proceed as follows: we fix p and, for the 
Gaussian profile, adjust ^jlq accordingly. A comparison of 
the solid lines in panels (a,b) of Figs. || and || shows that 
in the absence of transverse flow, at the same value of p, 
the dilution effects on (f)(p) from longitudinal expansion 
are almost identical for both density profiles. However, 
the additional flattening effects from transverse flow are 
weaker in the Gaussian model than for the box profile. 

We have demonstrated that the thermal model of ultra- 
relativistic heavy ion collisions does not predict an aver- 
age phase-space density characterised by a Bose-Einstein 
distribution function. Comparing data for ex- 
tracted from Eq. (^|) with such a distribution is likely 
to severely underestimate the pion chemical potential 
or temperature. The previously used ad hoc introduc- 
tion of a naive transverse blue-shift factor P,f3|, with- 
out first properly accounting for the dominant effects 
from longitudinal expansion, is based on faulty intuition 
and misleading. Longitudinal expansion strongly reduces 
the average freeze-out phase-space density at low p±, 
whereas transverse flow leads to an additional flattening 
of (/) (p±). At lowp^, where (f)(p) and thus the danger 
of distortions due to strong multi-boson symmetrization 
effects are largest, the strong longitudinal expansion 
of the heavy-ion collision fireball dominates the suppres- 
sion of the average phase-space density. An extraction of 
the pion chemical potential should be based on a compar- 
ison of the measured values for (f)(p) with Eqs. ([l2]), or a 
suitable generalisation thereof which more accurately ac- 
counts for shortlived resonance decays (which we treated 
here rather superficially). 
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